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Abstract 

Based on tests of a tentative detector for observing geo-synchrotron hard X-rays generated by primary electrons, we 
study the feasibility of probing cosmic electrons above a few TeV to over 10 TeV. Such high energy electrons are 
expected to give proof of sources near the Earth (e.g. supernova remnants such as Vela: age < 10 5 years located 
within <lkpc). The idea itself is rather old; a high energy electron emits synchrotron X-rays successively in the 
geomagnetic field and thus gives several X-rays aligned on a meter scale. This feature is a clue to overcome the 
background problem encountered in other traditional observation methods. We critically examine the feasibility of 
this approach assuming a satellite altitudes observation, and find that it is difficult to derive a precise energy spectrum 
of electrons but is possible to get a clear signal of the existence of several TeV electrons if the flux is comparable to 
the level predicted by a class of plausible models. For such observations, an exposure > lm 2 year would be needed. 
It would be attractive to incorporate the present scheme in the gamma-ray burst observations. 
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1. Introduction 

High energy electrons in space suffer synchrotron and 
inverse Compton scattering energy loss proportional to 
the square of their energy; this means high energy elec- 
trons cannot come from distant sources. If their energy 
exceeds several TeV, they are expected to be accelerated 
in and come from supernova remnants (SNR) of age less 
than ~ 10 s years, located within 1 kpc from the Earth. 
Among the small number of known SNR's satisfying 
such criteria, Vela is a good candidate yj]. Detection of 
such high energy electrons will give strong support for 



"By electrons, we mean electrons + positrons in this paper. 
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particle acceleration in near-by SNR and will give us 
important additional information related to high energy 
gamma ray observations, despite the fact that the elec- 
tron arrival direction cannot tell the source direction due 
to deflection by magnetic fields. 

Observation of high energy electrons is more difficult 
than many other primary cosmic rays due to the small 
flux and large background; we need large and thick de- 
tectors to cope with small flux and high energy, and spe- 
cial features in detector design to overcome the huge 
background mainly produced by protons. Therefore, 
experimental data over ~ TeV is very limited. Figure 
1 summarizea recent observations and some model pre- 
dictions. 

Recently, the HESS collaboration reported results in 
the TeV region by their air Cherenkov light detection 
method fl^yD. This method is quite different from other 
conventional direct observation methods and overcomes 
the small flux problem. The HESS result indicates a 
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steep spectrum, but it is probably premature to conclude 
the non existence of several TeV electrons^. We will 
need further investigation. 

The CALET projectJl and AMS-02 project|H are 
being prepared for observation on board ISS (Interna- 
tional Space Station). CALET uses a thick calorimeter 
and will be able to give a conclusive result beyond TeV 
after a few years of observation. However, it may turn 
out that we need to clarify the 10 TeV region flux (for 
which CALET may still be too small) to discuss near-by 
sources. AMS-02 uses a magnet spectrometer and may 
be limited for the TeV scale electrons. 

The idea of observing geo-synchrotron X-rays emit- 
ted successively by a TeV region electron and appear- 
ing almost on a line, for investigating cosmic primary 
electrons, is not new[8fl We study the feasibility of 
this idea, paying more attention to the background and 
stochastic nature of a realistic observation than earlier 
works. 

The CREST project [9] which is based on the same 
idea is preparing for a balloon height observation 
around Antarctica. Its result is highly awaited, and we 
will be able to see consistency with, or difference from, 
our estimation in the present paper, although we will 
discuss a satellite height observation. The CREST ob- 
servation may not be long enough due to the balloon 
flight limitation; nevertheless we can expect it will dis- 
close many new features in the geo-synchrotron obser- 
vation. 



2. Characteristics of geo-synchrotron X-rays from 
TeV region electrons 

Figure [2] shows a schematic view of geo-synchrotron 
X-rayfl We briefly recapitulate the characteristic fea- 
tures of geo-synchrotron X-rays which are produced by 
> TeV electrons. As is well known, the emissivity of 
synchrotron radiation by an electron of energy E e shows 
a peak at ~ y 2 v c where y = E e /m e is the gamma factor 
of the electron with mass m e , and v c the cyclotron fre- 
quency of the electron (in h = c = 1 unit). At higher 
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"It is interesting to note that their highest energy points are rather 
consistent with model A of the Vela source (Fig^ 

3 J. Nishimura had such an idea around the same time as this ref- 
erence and recommended us to explore the possibility. Although not 
easily accessible, there are even earlier works J6J (7J as mentioned in 
®. 

4 We use Cosmos (generation of X-rays)/Epics(detector response) 
for the present M.C simulation: http://cosmos.n. kanagawa-u.ac.jp 
The synchrotron X-ray emission is managed by Baring's formula 1 19] 
which can be used for our conditions as well as extreme conditions 
(very high magnetic field and very high energy). 
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Figure 1 : High energy primary electron spectrum obtained by Fermi- 
LAT: Aticl-2:[nll, AMS-2001:[Hl, ECCiO, HESS-2008:|1, 
HESS-2009:(l 1 BETS:[ll, PPB-BETS:OjJ. Model calculations la- 
beled galactic:| 18] and galprop:| 17] are for the sum of distant contri- 
butions. Others:|l] are possible contributions from near-by sources 
(A, B and C are used in the present paper). 



X-ray energies, there is a sharp exponential cutoff in 
the spectrum. The peak emissivity for the Earth reaches 
~10 keV for electrons of energy ~1 TeV. Therefore, if 
we observe several X-rays of energy 10 keV or more, we 
may expect that the average energy of X-rays, <Ex>, is 
related to the peak energy. If the average is proportional 
to the peak energy, E e cc m e y/<Ex> /v c is expected, but 
this seems to hold only when the minimum X-ray en- 
ergy is much smaller than the peak energy. 




Figure 2: Schematic view of geo-synchrotron X-rays. Enlarged in- 
layed figure shows an electron (leftmost) and accompanying X-rays 
with the density of several/m. 



Figure[3] shows the average integral energy spectrum 
of X-rays emitted by 1~20 TeV electrons which start 
isotropically from a height of 3 x 1 4 km and are directed 
to height 390 km, latitude 0°, longitude BO'S 



5 We used a geomagnetic filed of year 2010 by referring to the year 
2005 IGRF data; http://www.ngdc.noaa.gov/IAGA/vmod/igrf.html 
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Figure 3: Integral energy spectrum of geo-synchrotron X-rays emitted 
by 1 TeV to 20 TeV electrons 
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Figure 5: Integral production height distribution of those X-rays that 
have another X-ray within 50 cm. 
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Figure 4: Average energy of X-rays vs electron energy. The minimum 
X-ray energy is 1, 20, 100 keV as indicated in the Fig. For minimum 
energy of 20 keV, results for latitude 0° and 50° are shown. 



The number of X-rays and the spectrum shape de- 
pends on the location. At higher latitudes, in general, a 
harder spectrum is seen, and the total number of X-rays 
becomes smaller. However, such differences can be ne- 
glected for our present feasibility study, and we use this 
observation point for further discussions in this article^. 

In Fig. |H the average X-ray energy is plotted as a 
function of the electron energy. The minimum X-ray en- 
ergy is varied from 1 keV to 100 keV to see the change 
of the energy dependence. The result for latitude 50° is 
also shown to illustrate the small difference from the lat- 
itude 0° case. For 10 TeV electrons, the average X-ray 
energy is 600~700 keV. If the number of X-rays is lim- 
ited to several, the average generally becomes smaller. 



The background depends much more strongly on the latitude so 
we consider it separately. 



The integral production height of X-rays is shown in 
Fig. |5] We chose those X-rays that have another X- 
ray within 50 cm. The production height less than 500 
km contains 60 % of the events. To judge that X-rays 
observed are due to the geo-synchrotron emissions, we 
impose the condition that several X-rays come simulta- 
neously on a line. Since the emission angle of the syn- 
chrotron X-rays is order of ~ 1/y, the deviation from a 
line alignment is ~ 2.5 cm for X-rays coming from elec- 
trons of energy E e = 10 TeV at height 500 km. If we 
impose several X-rays to be aligned in a small distance, 
lower altitudes with higher magnetic field strength be- 
come more important. Therefore, 5 cm is a good mea- 
sure for the position resolution of the X-ray detector. 

In Fig. [6] we show the multiplicity distribution of 
X-rays falling in an area of 1.4 mxl.4 m for a primary 
electron of energy 10 TeV. The minimum energy of X- 
rays is set to be 40, 60 and 80 keV. In the distribution, 
we exclude the case where the electron itself enters the 
area; this is because the electron emits a number of pho- 
tons by radiation in the detector and they are observed 
simultaneously at various positions together with the X- 
rays. This makes it difficult to judge the intrinsic X-rays 
as discussed later (section [3~3l Fig. [7|l 

3. Consideration of the detector efficiency and back- 
grounds 

For further analysis, we have to assume a model de- 
tector and stochastic nature of the observation, together 
with a background for the observation. Our assumptions 
are based on a basic test experiment summarized in the 
Appendix. For example, we put an error to the energy, 
E (keV), absorbed in BGO by a random sampling with 
FWHM error of 16 V662keV/£ %. 
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Figure 6: Distribution of X-ray multiplicity in 1.4 m X 1.4 m area for 
a primary electron of 10 TeV 



3.1. The tentative detector 

The tentative detector is summarized as follows: 

1. The unit X-ray detector is made of BGO. The 
cross-section is 5 cm x 5 cm and the thickness 2 
cm (1.79 r.l; lr.l = 1.12 cm). Other BGO proper- 
ties are: density 7.13 g/cm 3 , light yield 8.2 pho- 
tons/keV or 15 ~ 20 % of Nal, decay time 300 ns. 

2. Each BGO is wrapped in teflon. The top part is 
covered by CFRP (Carbon Fibre Reinforced Plas- 
tic) of 500 jt/m thick. 

3. The gap between BGO units is filled with thin (1 
mm) Sn. Otherwise, a high energy albedo gamma 
ray could run the small gap rather a long distance 
and cause multi-Compton scattering which could 
be a major background. 

4. The whole detector size is 1 .43 m x 1 .43 m (28 x 
28 BGO's) 

5. Each BGO signal is read by a PMT at the bottom. 

6. Simple supporting platform is assumed together 
with material representing PMT (Fe and glass). 

3.2. The observation conditions 

We impose the following cuts on each event. 

1 . We put a threshold, E,h on the observed energy, Ex, 
in each unit BGO. The standard value of is 80 
keV. Those BGO's with observed energy Ex < E,h 
are neglected. 

2. The number of such BGO units (i.e, Ex > E t h), m, 
must be > 5. (Of course, we cannot know multiple 
incidence of X-rays in one BGO and such X-rays 
are regarded as a single X-ray, though the proba- 
bility is very small). 



3. The energy sum of the two highest energy BGO's 
(E\ + £2) must be < 4.5 MeV. This is to suppress 
background events due to multi-Compton scatter- 
ing of high energy gamma rays. 

4. Let the maximum distance among ra-BGO's be 
d cm. The center of each BGO is used to calcu- 
late the distance, d is converted to an integer b 
expressing an effective number of BGO blocks by 
(5.1/? + 5) = d. We require b > 11 (See section 
13.6.61 ). Hereafter b6 means b = 6. 

3.3. Collection power 

One of the merits of observing geo-synchrotron X- 
rays has been believed to be that we don't need to 
observe electrons directly and hence the effective area 
could be larger than that for direct observation. We de- 
fine the collection power, C p , by 



C„ 



Number of events with a desired condition 
Those events containing the electron 



-d) 



We note that we must exclude an event which con- 
tains the electron itself, since the high energy electron 
(> TeV) splashes a number of photons to surrounding 
BGO's and it becomes difficult to identify synchrotron 
X-rays (Fig. |7J. 
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Figure 7: The electron splashes a lot of photons to trigger a number 
of BGO's 



This fact seems to have been overlooked in past pa- 
pers, and unfortunately, reduces C p a large amount. If 
we use a loose cut for observation, C p could be as large 
as 5, but our condition reduces C p to < 1. This is partly 
due to the fact that we require m > 5; such high multi- 
plicity X-rays tend to be generated near the Earth where 
the magnetic field is strong, and hence separation of the 
X-rays and the electron is not large; that means many 
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of the high multiplicity events are accompanied by the 
electron. 

For our observation conditions, events containing the 
electron are automatically excluded since the energy 
loss in BGO's exceeds 18 MeV and E x + E 2 > 4.5MeV. 
In the present paper, we don't pursue how to utilize such 
events; if we could find a method for utilizing the events, 
C p will increase substantially. 

3.4. Dependence of observables on electron energy 

We use the average of square root of observed energy 
as a representative of m-BGO triggered by X-raysfl 
The energy is in keV unit. 



1.4mxl.4m 



V^/keV: 



(2) 



The distribution of e for various electron energies are 
shown in Fig. [8]for m > 5, £,/, = 80keV, and b\2 (see|4] 
in sec J3.2l > . 

The peak position, e p , of each e distribution and the 
electron energy are connected roughly by E e ~ (e p /lO) 3 
TeV. The width of the distribution is large (~ 80 % in 
FWHM) so that it would be difficult to derive a precise 
energy spectrum of electrons from those observations. 

However, there is a sharp exponential-like cutoff at 
the high energy side of the synchrotron X-ray spectrum 
and this gives us, in a sense, the function of a threshold 
counter to our detector. Actually, we see there is no 
sensitivity to electrons with energy lower than 2 TeV 
for the standard cuts. In Fig. [8] each distribution is 
normalized to give C p when it is integrated. The values 
are summarized in Table Q] 

Table 1 : Dependence of C„ on electron energy 



E e (TeV) 


3 


5 


7 


10 


15 


20 


c P 


0.0005 


0.018 


0.06 


0.13 


0.18 


0.18 



3.5. Case study of plausible models 

We use 3 models, A, B fljj] and C, for the electron 
energy spectrum in the highest energy region as shown 
in Fig. Q] 

Vela is a good candidate for high energy electron 
sources. Some models predict that ~ 10 TeV elec- 
trons are reaching the Earth. The prediction depends 



7 If we use a loose cut, the square root of X-ray energy is roughly 
proportional to the electron energy. However, this merit is lost by 
various cuts imposed on events. Using the simple average or median 
value of energy would lead to the same conclusion as the present one. 
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Figure 8: e distribution satisfying the observation conditions. Upper: 
for electron energy,£ c = 3, 5, 7 and 10 TeV. Lower: E e = 10, 15 and 
20 TeV. 



on the distance to Vela, diffusion constant etc. For 
example, distance dependence is shown in Fig. Q] for 
0.25, 0.30 and 0.35 kpc and a diffusion constant of 
D = 1.0 x 10 29 (£'/TeV) a3 cm 2 /s. An estimate of the 
probable distance to Vela is 0.25 kpc while the maxi- 
mum value is expected within 0.39 ±0.1 kpc [20]. Since 
other factors also affect the flux from Vela, we employ 
here the rather conservative (with respect to flux) dis- 
tance of 0.3 kpc and use the flux including the contribu- 
tion from Cygnus Loop. We call this model A. Larger 
diffusion constants flatten the spectrum. As an exam- 
ple of such a case, we use model B as shown in Fig. Q] 
Model C is a completely artificial spectrum to examine 
the detection limit. 

The line labeled 'galactic' is a model calculation lfl8ll 
which includes distant sources as well as possible near- 
by sources such as Monogem. The line labeled 'gal- 
prop' is another such model fl 1 711 calculated by the Gal- 
prop program. 

Table|2]shows the integral flux for model A, B and C. 
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Table 2: Integral flux of model A, B and C: electrons/;: sr-m -year 



Energy(TeV) 



model 


> 1 


> 3 


> 5 


> 7 


> 10 


> 15 


A 


2248 


243 


129 


74 


32.8 


8.1 


B 


3021 


324 


112 


50 


18 


3.6 


C 


2080 


86 


30 


14 


6 


1.6 



23D and EGRET [24] [25] data at higher energies (Fig. 



Chance coincidence by CXB must be considered at 
low energies (< 500 keV). As to low energy electrons 
and protons, their flux [ 27 ] [ 28 ] is much lower than CXB 
and we don't need to consider their effect. Also it's diffi- 
cult for high energy charged particles to generate BGO 
signals on a line with energy deposit less than a few 
MeV which is our maximum cut. 



3.6. Background considerations 

When we take an event with m(>5) synchrotron X- 
rays aligned, possible backgrounds could come from 

1 . Chance coincidence of X-rays or charged particles 

2. Multiple Compton scattering by a high energy X- 
or gamma-ray entering the detector 

Possible sources of such backgrounds would be 

1 . Uniform isotropic X-rays and gamma-rays ( Cos- 
mic X-ray Background or CXB) 

2. Cosmic ray charged particles such as electrons and 
protons 

3. X-rays from the galactic plane and strong point 
sources 

4. Albedo particles from the earth's atmosphere 

5. A Gamma ray burst 

6. Radiation from the environment surrounding the 
detector 

7. South Atlantic Anomaly Radiation 

8. A Solar flare 

As to the environmental radiation, we must choose a 
low background environment and, this is expected to be 
possible owing to our observation conditions. The last 
two sources might be avoided by switching off the elec- 
tronics. If a gamma ray burst is as strong as triggering 
more than 5 BGO's almost on a line in a short time (10 
ns; see later), we could get information from other GRB 
observatories. 

Next, we summarize other factors and derive the 
chance coincidence rate. 

3.6.1. Cosmic X-ray Background: CXB 

CXB (Diffuse X-rays) is now well understood and we 
can calculate the background due to CXB. The back- 
ground by other sources is estimated as the ratio to the 
one for CXB. As the low energy CXB, we use the data 
from HEAO-1 112 111 which was confirmed by Swift/BAT 
1I22I1 recently. The data can be continued to the SAS-2 



3.6.2. Atmospheric albedo particles 

High flux low energy particles which might con- 
tribute to chance coincidences can be absorbed by the 
support structure to be put on the Earth side. High en- 
ergy particles are not able to make fake events by multi- 
Compton scattering. However, X-rays and gamma rays 
in the intermediate energy region must be considered. 
Their angular distribution and latitude dependence are 
rather complex. At hi gh e nergies (> 30 MeV), there 
are data from SAS-2 Q and EGRET Q. In the 
X-ray energy region, recent reliable observations from 
Swift/BAT [22] show a steeper spectrum and thus a 
higher flux in the low energy region ( < 200 ke V) than 
the older one 113111 . 
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Figure 9: Energy spectra of the isotropic diffuse X/gamma-rays 
(CXB) and atmospheric albedo X/gamma-rays. 



Although there is a maximum of 5 times difference 
in the fluxes around the equator and magnetic poles, 
the over all difference is within a factor of 3. We 
smoothly connect Swift/BAT data averaged over lati- 
tude and SAS-2(~ 30 MeV) and EGRET (~ 100 MeV) 



8 Although energy is too high for our interest, a recent Fermi- LAT 
[2g] result over several 100 MeV region gives a lower intensity than 
12511 which gives little bit lower revised flux than the original one l24h 
at several 10 MeV. 
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Table 3: Flux of CXB: the 3rd line shows the integral flux over ^(in 7rsr) 



E x (MeV) 


0.01 


0.02 


0.03 


0.04 


0.05 


0.06 


0.07 


0.08 


0.1 


0.2 


/cm 2 s sr MeV 


316 


102 


42.2 


25.6 


15.0 


9.3 


6.15 


4.33 


2.43 


0.417 


/cm 2 s 


10.7 


5.18 


3.12 


2.09 


1.47 


1.10 


0.86 


0.704 


0.5 


0.175 



0.3 


0.5 


1 


2 


5 


10 


20 


50 


100 


0.15 


0.042 


7.4(-3) 


1.4(-3) 


1.5(-4) 


2.8(-5) 


5.9(-6) 


7.9(-7) 


1.6(-7) 


0.096 


0.045 


1.69(-2) 


6.6(-3) 


1.75(-3) 


7.2(4) 


3.08(4) 


1.01(4) 


4.5(-5) 



(-2) etc means 10 2 etc 



data. Figure [9] shows the energy spectra of CXB and 
albedo X/gamma-rays which we used in the simulation. 
At the high energy side, the albedo flux is much higher 
than CXB. Since the atmospheric albedo is complex, we 
use, for safety, an overestimated flux if there are uncer- 
tain factors. 

From this point of view, we employ a 1 .5 times higher 
flux than Swift/BAT data in Fig. [9] We also assume that 
the Swift/BAT data is the flux from the Nadir direction 
(this will probably lead to overestimation of the back- 
ground). If we measure the X/gamma-ray angle 6 from 
the Nadir, the Earth horizon is at 6 ~ 70° where the 
albedo flux becomes maximum. This enhancement ef- 
fect can be well modelled by assuming that the angular 
distribution of X/gamma-rays which enter the bottom 
side of the detector is not cos 6d cos 9 but dcosO (i.e, 
not isotropic but there is ~ sec enhancement). 

With this assumption, if we compare the all angle in- 
tensity of the albedo over 80 keV to CXB without con- 
sidering absorption by the structural material, we get 
albedo/CXB ~ 4. If we assume a supporting material 
equivalent to Pb 500 /im or Fe 1 mm, the albedo flux be- 
comes 1.29 (/cm 2 -s) or 2. 15 (/cm 2 -s), respectively. The 
ratio to CXB reduces to ~ 1.8 or ~ 3.1, respectively. 
However, we don't consider this absorption effect for 
albedo for the moment. 

3.6.3. X-rays from the galactic plane 

The X-ray intensity per unit solid angle from the 
galactic plane is factor ~10 stronger than CXB. How- 
ever, what is important for the chance coincidence is 
the integral value over a hemisphere. Its ratio to CXB is 
considered to be almost constant above 10 keV to over 
10 MeV. This ratio is estimated to be less than 0.1 by 
Ginga[32] and ASCAJH]. This value is for the case 
when the galactic plane is near the zenith. Therefore, we 
can safely neglect X-ray contribution from the galactic 



plane as compared to albedo or CXB. 

3.6.4. Chance coincidence 

We first consider the chance coincidence by CXB; we 
estimate background events with m X-rays aligned. The 
factors to be considered are 

• Time resolution of the system, r in FWHM. 

• Effective flux over 80 keV. / 

• Area where the X-rays there can be regarded as on 
a line. AS 

• The total area of the detector. S . 

• Observation time. T. 

• X-ray multiplicity, m. 

Then, the number of chance events, N, in the observa- 
tion time is 

N ~ fST(fAS2T) m - 1 (3) 

Expressing t in ns, area in m 2 , T in year and / = afo 
where /o = 0.704/cm 2 -s is the CXB integral flux over 80 
keV (we call a effective background index), we obtain 

N ~ 2.18aSr(aASry"- 1 (1.408xl0- 5 ) m ~ 1 xl0 11 (4) 

For S = 2, AS = 0.2, T = 1, m = 5, 

TV- 2.74 x 10 _11 ar 5 T 4 (5) 

This relation is shown in Fig. [10] 

As described later, we found r = 5.4 ns with use of 
a constant fraction discriminator; this value is realistic 
even at 80 keV but we assume here r =10 ns for safety. 
The maximum possible value of a would be 

a = l(CXB) + 4(albedo) + 0.1 (galactic plane) ~ 5 
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Figure 10: The number of chance coincidence events, N, as a function 
of time resolution (r) and effective background index(a) for S=2m 2 , 
T = 1 year 



Then, we can expect the number of chance coinci- 
dence events is <s 0.01 for a 1 year observation. As the 
condition for line alignment, we use only the AS con- 
straint, so the actual background could be reduced. 

3.6.5. Background by multi-Compton scattering 

A gamma-ray may repeat a Compton scattering in 
BGO and/or environmental media and can make a fake 
event as if synchrotron X-rays are on a line. Due to the 
condition, m > 5, the energy of such a gamma-ray is 
concentrated in 5 ± 3 MeV and very seldom exceeds 25 
MeV. 

If we require m > 5 and the maximum distance be- 
tween the triggered BGO's be greater than 30 cm, we 
may expect a very small number of fake events. How- 
ever, the number of CXB over 1 MeV entering 2 m 2 area 
in one year exceeds 10 10 while expected signal is order 
of 10 ~ 20 events. Then, we have to consider an event 
with probability of ~ 1/10 9 . 

We first calculate multi-Compton event rate using 
CXB, and the contribution from albedo is considered 

. . , . . _ . . , albedo 

by introducing the effective ratio a = 

3 5 CXB 

Multi-Compton events are efficiently produced by ~ 

5 MeV gamma-rays. The albedo gamma-rays in this 
energy range are ~ 40 times more abundant than CXB. 
Albedo from the Earth horizon have a nadir angle of 

6 ~ 70 degrees and are apt to produce multi-Compton 
scattering events. Therefor, it would be appropriate to 
take a' ~ 40 x 3tt/4 ~ 100 for safet>0 



9 sin 8 is a good enhancement factor for multi-Compton scattering 
(8 is a nadir or zenith angle). Then, albedo enhancement factor rel- 
ative to CXB is J sin 8d cos 8/ j sin 8 cos 8d cos 8 = 3^/4 where the 
integration region from to 2n may be used. 




1)10 

bll 
1)12 



Model A 
19-20 events 



1.4mxl.4m 

m > 5 
Eth=80 keV 




Figure 11: e distribution by model A and background for one year 
observation. The background rapidly reduces as b increases. Top: £>6, 
background dominates, next: b&,b9, bottom: blO, 11, 12. background 
vanishes but signal remains almost the same. 
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3.6.6. Signal vs background 

To show how the aligned X-rays are widely spanned, 
we use, for example, MO as defined in[4]in section [3721 

Figure QT| shows signal vs background for model A 
under the condition of S = 2m 2 , E,h =80 keV, m > 5 
for various The CXB background for b9 becomes 
~ 0.2 but if the albedo effect is to be included, a' must 
be multiplied and the background becomes quite com- 
parable to the signal (~ 20 events). The signals for blO 
to b 12 do not change appreciably, while background de- 
creases faster than exponential and for b 10,1 1,12 it be- 
comes -0.02, 0.001, <K O.OOlB Therefore, with bll 
or bl2, we can completely neglect the background from 
albedo multi-Compton scattering. 



fe; 

T3 



A ~19 events 



e(=< VEx/keV) >) 
Figure 12: Comparison of expected signals by Model A, B and C 

Figure [12] compares the case of Model B and C with 
Model A (bl2). The number of expected events for 
Model B is ~ 14; it's smaller than model A (~ 19) 
but within the detectable range. The smaller number re- 
flects the difference of the spectrum shape, but it would 
be difficult to derive the original spectrum shape. Arti- 
ficial model C gives ~ 3.2 events in a year. This gives 
a reference of the minimum detectable flux; we would 
be able to say that a model is plausible or incompatible 
with the observation only if the flux is much higher than 
1.5 x 10- 8 (£/2.3TeVr 3 (/m 2 s sr GeV) over 2 TeV. 



'"Simulations corresponding to several years are converted to one 
year equivalent. In this simulation, only CXB with energy greater than 
1 MeV is considered. 

"This faster-than-exponential-decrease is realized only if we fill 
space between BGO's by, say, Sn. If there is 2 mm unfilled space be- 
tween BGO's, the event rate (/year) becomes ~ 300exp(-1.15(fo - 5)) 
and even bl2 will receive albedo background of order ~ 10. With 
the filled space, background rate for b > 6 is expressed by ~ 
300exp(-1.15(;t - 5)' 30 ). This relation is inferred by extrapolating 
the calculations done up to blO (calculations for bll and bl2 need 
huge compter power). 



3.6. 7. Light/thermal shielding and pile-up effect 

In this trial calculation, we assume X-rays with more 
than 80 keV energy. However, we must prevent a pile- 
up effect in BGO due to lower energy X-rays. CFRP is 
supposed for light/thermal shielding, and if we make its 
thickness to be 500/mi (0.11g/cm 2 ) carbon equivalent, 
we obtain the X-ray transmission rate as shown in Table 


CXB over 3 keV is expected in a 5 cm x 5 cm area at 
~400 Hz, but those photons reaching the BGO are less 
than 100 Hz. Low energy albedo has much smaller flux 
than CXB, as in Fig. [9] Supposing the effective flux 
of albedo over 80 keV is 4 times CXB, then the rate of 
each BGO is ~ 70 Hz. Therefore, X-rays entering each 
BGO will not exceed 200 Hz in normal conditions. This 
means pile-up in BGO dose not matter since the decay 
time of BGO light emission is ~ 300 ns. 

3.6.8. Trigger rate and dead time 

Trigger rate will be governed by multi-Compton scat- 
tering. According to the simulation, the CXB contri- 
bution with the condition of any 3 BGO coincidence 
(Ex > 80keV) is ~lHz, and for any 4 coincidence 
~0.045Hz. The albedo will contribute a maximum of 
100 times of this and will make ~100Hz for any 3 and 
~ 4.5 Hz for any 4 coincidence, respectively. Our tar- 
get is m > 5, however, any 4 coincidence trigger will be 
appropriate. 

Since the decay time of BGO light is ~ 300 ns, we 
may assume ~ 1ms dead time. Then, the trigger rate of 
~10 Hz presents no problem. 

4. Summary and concluding remarks 

• The purpose of this article is a study of the feasi- 
bility of probing cosmic ray primary electrons in 
the ~10 TeV region by observing geo-synchrotron 
X-rays. According to some plausible models, such 
high energy electrons are expected to come from 
supernova remnants such as Vela. 

• For this, we assumed a tentative detector design 
consisting of a number of BGO blocks as described 
inlXTl 



It is found that we cannot use events which contain 
the electron itself falling on the detector; a high 
energy electron splashes a number of X/gamma- 
rays and thus it is difficult to identify aligned X- 
rays. This reduces event the rate substantially. 

To realize a background free observation, the num- 
ber of triggered BGO's (with X-ray energy > 80 
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Table 4: X-ray absorption by 500/jm carbon 



E x (keV) 


3 


5 


10 


60 


100 


incident angle (deg) 





70 





70 





70 





70 





70 


transmission rate(%) 


0.0012 





25 


1.5 


75 


43 


98 


95 


98 


95 



keV) must be m > 5 and they must span more than 
62 cm. 

• It will be difficult to get a spectrum shape for elec- 
trons. However, owing to the exponential sharp 
cut-off of the X-ray energy spectrum for each elec- 
tron energy, the assumed detector has no sensitiv- 
ity to electrons below 2 TeV, and thus works like a 
threshold detector. 

• We can verify the plausibility of models if the 
flux is much (4~5 times) higher than 1.5 x 
10- 8 (£/2.3TeVr 3 (/m 2 s sr GeV) over 2 TeV. The 
number of signal events is expected to be 15 to 20 
in a year. 

• It will be possible to increase the number of signal 
events by lowering the threshold for observation 
down to sim50 keV without a big problem (chance 
coincidence is still negligible and multi-Compton 
scattering event rate is not affected). Then, the 
number of signal events will double or the detec- 
tor could be made smaller. 

• Lowering m by 1 is very challenging because the 
chance coincidence probability increase substan- 
tially. However, in our calculation we only consid- 
ered AS for a line alignment constraint. This is a 
loose condition rather than a rigorous one. Also we 
may note that e by chance coincidence distributes 
sharply at the lower end of the true signal distri- 
bution. This is useful to separate some of the true 
signal. These factors are worth further considera- 
tion. 

• The observation considered here would be more at- 
tractive if it could be combined with gamma-ray 
burst observations. 
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Appendix A. Basic test of detector components 

To get basic performance of assumed detector com- 
ponents, we performed the following basic test experi- 
ments. The parameters thus obtained were used in the 
present calculations. 

Appendix A.l. Materials 

• BGO (SICCASS 

5cmx5cmx2cm t block x8. As reflector, white 
teflon was used. Thin CFRP was used for 
light/thermal shield. 

• PMT (Hamamatsu): 

R3318-HA(bialkali, 2"square) x2. 
R6231-100HA(superbialkali, 2'» x2 

To see possible individual differences, two PMT's 
were used for each type. 

• X and gamma-ray source 

241 Am(60keV, 14 keV), 57 Co(122keV), 137 Cs(662 
keV), 60 Co(1.17MeV, 1.33MeV) 

Appendix A.2. Temperature dependence 

To avoid a temperature dependence in the test, and 
for future applications, we examined the temperature 
dependence of pulse height for the complete absorption 
peak of 662 keV gamma-rays. 

At 15 to 30° C, temperature dependence of BGO is 
within ~ +2%. At lower temperatures down to -30°C, 
0.5%/deg dependence is seen. Therefore, the tem- 
perature effect is negligible in our test conducted at 
room temperatures. LaBr3 has only ~ ±2% changes 
over the entire range. 



A LaBr3 (Sangoban) package was partly used as a reference, es- 
pecially, at low energies where BGO output becomes weak. LaBr3's 
property: density 5.29 g/cm 3 , radiation length 2.1, decay time 25 ns, 
light yield 63 photons/keV or 200 times Nal. 



10 



Appendix A.3. Energy resolution 



Appendix A.4. Dynamic range 



We exposed isotope sources mentioned in 



Appendix A.l to the BGO's to obtain pulse height 
distributions. Except for 60 keV X-rays from 241 Am, 
data was taken for all combinations of 3 sources, 8 
BGO blocks and 4 PMT's (total 96 cases). 

The maximum difference of energy resolution (in 
FWHM) among 8 BGO's is 4% for 662 keV and 7% 
for 122 keV. So we need to choose good BGO for low 
energy X-ray observation in actual application. Two 
PMT's, Bialkali (B) or Super Bialkali (SB), show only 
a small difference (within 2%). SB gives ~ 2 % better 
resolution than B for 662 keV and ~ 4% for 122 keV. 
The SB PMT recognized clearly two peaks at 1.17 MeV 
and 1.33 MeV gamma-rays from 60 Co, but the B PMT 
showed no clear separation. 

This difference is reasonable considering the maxi- 
mum yield wavelength (480 nm) of BGO, wave length 
dependence of the quantum efficiency and differences in 
light collection efficiency due to the square and circular 
PMT areas. 







1 1 1 1 
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10 100 1000 

Energy (keV) 

Figure A. 13: Energy dependence of energy resolution (% of FWHM). 
The results for SB PMT voltage of 900 V and 1000 V are shown. 
Accuracy of points with () is less than other points due to difficult 
background subtraction. The point at ~ 33 keV is estimated from a 
peak obtained in 57 Co case (Note: the K-shell energy of Bi is 90.5 keV 
and 57 Co has ~1 1% 136 keV line besides its 122keVline. We expect 
a line near (123.5 - 90.5) keV. ). The curve shows the resolution, 
16 V662 (keV)/E) %, used in the simulation to put a random error on 
the absorbed energy in BGO. 



Figure |A 1 3 1 shows SB PMT measurement results. A 
difference of ~3 % does not matter for our purpose and, 
in the simulation, we assumed a resolution shown by 
the curve in the figure to introduce a random error to the 
absorbed energy in BGO. 



Dynamic range of the PMT was examined by apply- 
ing voltages from 850 to 1000 V. The pulse height is lin- 
ear as a function of X-ray energy from 14 keV to 1.25 
MeV, except for 850 V case for which linearity is lost 
in the 10 keV region (it is difficult to see the 14 keV 
peak). Voltage dependence of the pulse height is oc y 5 6 
for each energy. Since our target is 50 keV to few MeV, 
this feature will make it easy to build electronics. 

Appendix A.5. Time resolution 



BGO 



X 



J PMT 




AMP 




DISC 












PMT 




AMP 




DISC 



stop 



BGO 



Figure A. 14: Setup for measuring time resolution of the BGO+PMT 

system 



To avoid chance coincidence, good time resolution 
is indispensable. To see the time resolution, we use a 
system as shown in Fig. IA. 141 and measured time dif- 
ference, AT, from start to stop. Two gamma rays from 
60 Co (1.17 and 1.33 MeV) are emitted isotropically and 
simultaneously and will enter BGO at the same time. In 
some cases, one gamma-ray which enters a BGO block 
may be Compton scattered and enters another BGO. 
This timing is also the same within our accuracy. 




1000 1200 
AT (channel) 



Figure A. 15: Time difference distribution observed in the setup shown 



in Fig. IA. 141 One channel is 50 ps. The constant fraction method 
gives 5.4 ns FWHM resolution, while the leading edge method gives 
9.1 ns, and the distribution has a long tail deviating from Gaussian. 



To get good resolution, we use a constant fraction dis- 
criminator (ORTEC 935). We also use a leading edge 
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discriminator to confirm the merit of using the constant 
fraction method as seen in Fig. IA. 151 

The pulse height from 1.2 MeV absorption is ~ 100 
mV, and we use a discrimination level of 50 mV. For 
timing, 50% fraction is used. The data include Compton 
events which have smaller energy than 1.2 MeV; thus 
the result indicates similar good time resolution even for 
smaller energy deposit. To see this, we use the discrim- 
ination level of 25 mV and 10 mV to find less than 1 ns 
degradation of the time resolution. This suggests 10 ns 
resolution is quite safe for several tens keV region. 
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